Angular resolved specific heat in iron-based superconductors: the case for nodeless 

extended s-wave gap 
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We consider the variation of the field-induced component of the specific heat C(H) with the 
direction of the applied field in _Fe— pnictides within quasi-classical Doppler-shift approximation, 
with special emphasis to recent experiments on FeSeo.4Teo.6 [Zheng et al., arXiv: 1004.2236 . We 
show that for extended s— wave gap with no nodes, C(H) has cos 40 component, where (j) is the 
angle between H and the direction between hole and electron Fermi surfaces. The maxima of C(H) 
are at ■7r/4, 37r/4, etc. if the applied field is smaller than Hq < IT, and at = 0,n/2, etc. if the 
applied field is larger than Ho- The angle-dependence of C(H), the positions of the maxima, and 
the relative magnitude of the oscillating component are consistent with the experiments performed 
in the field of 9T >> Ho- We show that the observed cos4(/!> variation does not hold if the s— wave 
gap has accidental nodes along the two electron Fermi surfaces. 

PACS numbers: 74.70.Xa, 75.10.Lp, 75.30.Fv 



The symmetry and the structure of the superconduct- 
ing gap in Fe— based superconductors have been subjects 
of numerous experimental and theoretical papers in re- 
cent yearsir— . There is a growing consensus among re- 
searchers that (i) the gap has an extended s— wave sym- 
metry - it belongs to a symmetric Aig representation of 
the symmetry group of a square lattice and its aver- 
age values along electron and hole Fermi surfaces (FS) are 
of opposite signs; (ii) that superconductivity originates 
from umklapp process in which pairs of fcrmions hop be- 
tween conduction and valence bands; and (iii) that the 
pair-hopping interaction is a combination of a screened 
Coulomb interaction and a magnetic interaction, medi- 
ated by antiferromagnetic spin fluctuations. 

A more subtle and currently hotly debated issue is 
whether the gap has nodes. This is not a symmetry is- 
sue as, quite generally, extended s— wave gap can be ap- 
proximated by a constant only along the hole FS, while 
along the two electron FSs it has angle-independent and 
008 2^/! components: Ae(?/;) = Ae(l + bcos2ip), where ip 
is the angle counted from the line connecting the two 
FSs, and the sign of Ae is opposite to the sign of the gap 
along hole FSs.- Such Ae{ip) has no nodes if & < 1 and 
has "accidental" nodes when & > 1 at non-symmetry se- 
lected directions cos 2^ = 1/6 Because i^e— pnictides are 
multi-orbital systems, the cos 2-0 component of the inter- 
action is generally not small, i.e., b can be either larger or 
smaller than 1, depending on the material. Furthermore, 
b gets larger when one includes into the gap equation 
intra-band Coulomb repulsion because this term couples 
to the gap averaged over the FS and hence reduces angle- 
independent gap components but does not affect cos 2-0 
components (jje). As a consequence, b becomes progres- 
sively larger as the system moves further away from the 
SDW phase and the effect of intra-band repulsion grows, 
that is, overdoped ferropnictides are more likely to have 
nodes in the gap. 

The issue whether or not the gap in i^e-pnictides has 



nodes is crucial for the understanding of low-energy prop- 
erties of these materials and deserves a careful study. The 
subject of this work is the interpretation of recent high- 
accuracy measurements^ of the dependence of the low- 
temperature specific heat in FeSeo.4Teo.6 on the direction 
of a magnetic field. Similar experiments have been car- 
ried out in the past on borocarbides^, and heavy-fermion 
CeCoIns (Ref. i) and revealed modulations generally 
consistent with the d-wave gap (for details seeiSiil; for 
experiments on thermal conductivity see^-). 

The generic reason for field-induced modulations of 
specific heat and thermal conductivity in unconventional 
superconductors is that a magnetic field induces vortices 
along the field direction. In a vortex state of a type-II 
superconductor, scattering of quasiparticles on vortices 
gives rise to a non-zero density of states (DOS) at zero 
energy. The magnitude of this residual DOS depends on 
the angle the field makes with the position of the minima 
of the modulus of the superconducting gap. This leads 
to modulations of the field-induced linear-in-T term in 
the specific hea t^°i^^'^'^ . This reasoning works the best 
when the gap has nodes but should be generally applica- 
ble also to materials where the gap varies along the FS 
but not necessary has nodes, provided that the field is not 
too small. Fe-pnictides are strong type-II superconduc- 
tors (both magnetic and coherence lengths are of order 
2 — 3 X 10^ A, much smaller than the penetration depth 
A ^ 3 X lO^A^'*), and vortex state extends to almost all 
fields [the upper critical field is about lOOT, lower critical 
field is about iJ^i < lOmT^^]. 

In FeSeo.4Teo.6, the datai show cos 40 modulation of 
C(H), with the magnitude of about 1% of the total field- 
induced C(H). The maximum of C(H) is at = 0,7r/2, 
etc. what correspond to the directions of H along the 
axis between hole and one of electron FSs in the un- 
folded Brillouine zone (BZ) [along the diagonals in the 
folded BZ]. The cos 40 modulation of C(H) was origi- 
nally interpreted^ as evidence for the nodes in the gap. 



However, to be consistent with the observed near-perfect 
cos 40 form of C(H), the nodes have to be located pre- 
cisely at 45° with respect to the x axis, i.e., right at the 
crossing points of two electron FSs in the folded BZ (see 
Fig. [1]). This is generally incompatible with the "acci- 
dental" nodes located at some arbitrary angles 0. The 
authors of Ref. argued that the data are inconsistent 
with no-nodal extended s— wave gap and an extended 
s— wave gap with accidental nodes. To explain the data, 
they included spin-orbit coupling and argued that it cre- 
ates nodes on electron FSs at exactly 45°, even if the gap 
was nodeless in the absence of spin-orbit interactionii. 

In this communication, we argue that the data of 
Ref. can be actually explained quantitatively in a con- 
ventional semi-classical Doppler-shift scenario for field- 
induced oscillations of C (H) , but only if one assumes that 
the gaps along the electron FSs have no nodes. The node- 
less gap in FeSeo.4Teo.6 has been extracted from STM 
data^^, and we argue that STM and angle-resolved spe- 
cific heat measurements are consistent with each other. 

Our reasoning is two- fold: First, as we said, the two 
electronic gaps generally have the forms Ae = A(l ± 
6 cos 2-0). The formula for the specific healji^ contains 
Ag, i.e cos 20 and cos 40 terms. The cos 20 terms can- 
cel out when the contributions from the two electron FSs 
are added, while the cos 40 term generates cos 40 mod- 
ulation of C(H). This, however, holds only if Ag does 
not have nodes {h < 1), otherwise the modulation of |Ae| 
will be more complex leading to a more complex struc- 
ture of C(H). Second, to be consistent with the data, 
the sign of cos 40 modulation of C (H) should be positive 
(maxima should be at = 0,7r/2, etc). At small H this 
is not the case - the sign is negative. We show, how- 
ever, that the sign of the cos40 term in C(H) depends 
on the magnitude of the field and changes from negative 
to positive as the field increases. We estimated the field 
where the sign changes and found that it is about IT for 
all & < 1, much smaller than 9T, at which experiments 
have been performed. In other words, at the field of 9r, 
cos 40 oscillations of C(H) have maxima at = 0, 7r/2, 
etc for arbitrary strong oscillating component of Ag, as 
long as it remains nodeless (i.e., as long as 6 < 1). Fur- 
thermore, at 9T field, the magnitude of the oscillating 
part of C(H) is around 1% of the total specific heat, like 
in the data^, and this number weakly depends on h ex- 
cept for very small values, where it vanishes as . The 
conclusion of our analysis is that the data on C(H) are 
quite consistent with the "conventional" theory of field- 
induced oscillations, provided that the gaps along elec- 
tron FSs have moderate cos 20; oscillations and no nodes. 
The range 6 « 1 is a gray area, and in the presence of 
some amount of disorder modulations of C(H) may still 
look like cos 40 even when the gap has pairs of weakly 
spaced "accidental" nodes. Still, a more natural expla- 
nation of the data in FeSe^ j^Te^.d is that the gap has no 
nodes. 

The sign change of the prefactor for the oscillating 
cos 40 component in the specific heat and thermal con- 
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FIG. 1: (color online) Fermi surface topology of iron-based 
superconductors in the unfolded (a) and folded (b) BZ. There 
are two small hole pockets centered around the F— point and 
two elliptic electron pockets centered around the (ivr, 0)- 
point and (0, ±7r)-point of the unfolded BZ. Upon folding two 
electron pockets fold into the (tt, 7r)-point of the folded BZ 
(dashed-dotted line in (a) is the boudary of the folded BZ). 
The magnetic field rotation is shown by angle 00- 

ductivity is the well-known phenomenon for d— wave su- 
perconductors. The detailed theoretical study of the sign 
variation of the prefactor for the cos 40 term with chang- 
ing magnetic field and temperature has been performed 
by Vorontsov and Vekhter (VV)^^ and by Hiragi et alii. 
VV recently performed numerical analysis of the angular 
dependence of C(H) in the ferropnictidesi^, and found the 
change from negative to positive prefactor of the cos 40 
term with increasing field and temperature. Our results 
are fully consistent with theirs and provide analytical rea- 
soning for the sign change of the cos 40 term in the iron- 
based superconductors. 

Like we said, scattering of quasiparticles on vortices 
gives rise to a finite DOS at zero energy what in turn 
gives rise to a linear in T specific heat: 

C(H) - AT7V(H), N{li)= j^'j dioY^N^i^A 

(1) 

where A is the overall factor and j is the band index. The 
experiment^ has been performed at the low T ^ 2.6K 
when terms of higher order in T are irrelevant. 

We consider FS geometry consisting of two hole FSs 
centered at F point and two electron FSs centered at 
(0,7r) and (tt, 0) in the unfolded BZ, or at {tt,tt) in the 
folded BZ (Fig. [T|). The potential presence of the third 
hole FS is not essential for our analysis because oscil- 
lations of C(H) come only from the two electron bands. 
For simplicity, we assume that all bands are circular, i.e., 
neglect ellipticity of electron bands. We will follow Refji^ 
and employ the formula for N^{Ii,^j) obtained by solv- 
ing semi-classical Eilenberger equations for a given vor- 
tex lattice within Brandt-Pesch-Tewordt (BPT) approx- 
imation in which the dependence on the normal Green's 
function on the center of mass coordinate of a pair is re- 
placed by an average over a unit cell of the vortex lattice. 
Hiragi et alii computed the DOS beyond BPT approxi- 
mation and obtained, but the changes turned out to be 
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FIG. 2: Calculated angular dependence of the density of 
states for large a — 10(a) and small a = 0.1(b). We set 
b = 0.5. The solid and dashed curves refer to the numerical 
solution of Eq.(3) and the approximate analytical formulas 
Eqs. (4) and (5), respectively. In panel (a) we matched ana- 
lytical and numerical curves at tf) — 0. 
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(color online) (a) The functional form of the function 
from Eqs. (6) and (7). Sign change of F{a) implies 
sign change of the prefactor for cos4(^ term in the density 
of states; (b) the difference between A'^(<^) for (f> = 7r/2 and 
(j> = -K /A from Eq. (3) as a function of a for various h. For 
any 6 from the interval < 6 < 1, N{n/2) — N{tt/A) changes 
sign at a finite a = ao, ranging between 0.65 and 2.8. 



quite small. Suppose that H is applied in Fe — Fe plane, 
at an angle (f> with respect to the x— axis (which in mo- 
mentum space is the direction between hole and electron 
bands). In the area surrounding the vortex, the DOS at 
zero energy can be generally written as^ i^^'^'^ 

iV-(H,^)^ ^ (2) 
y'a2(H,v.) + (A^-(V'))' 

where a(H, ip) is proportional to the component of the 
Fermi velocity normal to the field Vp — vpsm{ip — (f)): 
a(H, 0) = asin{'i(j — 0), where a = cwf/(2\/2A), A — 
yjhc/ (2|e|B) is the magnetic length, and c — 0(1) is 
a numerical factor which carries information about the 
geometry of the vortex structure and the distance from 
the vortex core^. 

Substituting Ag = A(l ± 6 cos 2-0) and A/i = const 
into ([2]) and shifting the integration variable, we obtain 
iV(H) = Nh{H) + A^e(H), where Nh is independent on 
the direction of H, and A^e(H) g(iVe(H, V) + 

A^e(H,V' + 7r/2), where 



V ^ ^ 

and = 4^ = Hq/H. 

Consider the limits of small and large a separately. At 
large a (small fields) A^e(H) cx 1/a. Oscillating com- 
ponent of Nf. can be obtained analytically at small b. 
Expanding in b we find 

N,{ll) = N,{H)--—-cos^ + 0{b^) (4) 
157r|a| ^ ' 

where Ne{H) is a 6— dependent non-oscillating term. We 
see that the DOS does contain cos4(/) oscillations, as we 
anticipated, however the sign of the cos Acj) term is nega- 
tive, which implies that C(H) has peaks at 7r/4, 37r/4, 
etc, in disagreement with the data. At larger b, the 
oscillating part of A^e(H) contains higher harmonics 



cos 80, cos 12(/)... and has to be calculated numerically. 
We present the results in Fig. [DJa). We see that the 
oscillating part of iVe(H) still well described by cos 40 
form for arbitrary 6 < 1, despite that higher harmonics 
are not parametrically small. Also, the sign of the oscil- 
lating part remans negative at large a for arbitrary & < 1 
(i.e., for all 6 < 1, the maxima of C(H) are at 7r/4, etc.). 

The situation changes, however, in the opposite limit of 
large fields, when a << 1. Now iVe(H) can be expanded 
in a. The expansion requires care because of infra-red 
divergencies and yields, at arbitrary b < 1 

A^e(H) ^2- — + 2a262cos40 (l + 0(a)) + ... (5) 

TT 

where dots stand for terms of order a^, at which order 
higher harmonics appear. We see that the oscillating 
component is now cos 40 for all 6 < 1, and the sign of the 
oscillating part is positive, i.e., the maxima of C(H) are 
now at = 0,7r/2..., like in the experimental data^. In 
Fig. we present the result of numerical evaluation of 
Ne{il) and compare it with Eq. ([5]). Clearly, there are 
cosAcf) oscillations with a positive prefactor. The value 
of a at which the oscillating part of A^e(H) changes sign, 
and the crossover from a small field to a high field behav- 
ior can be analyzed analytically at small b. Expanding 
Eq.(3) in b to order b^ and integrating over ip we obtain 

7Ve(H) = - arctan - + a^b^ cos40i^(a) (6) 
TT a 

where 

^/ N 1 , > sini/'cos4-0 /„ 9 . 9 ,N 
F{a) = — / di> , ; {2a' - sin^ V (7) 

We plot a'F{a) in Fig. ISja). This function changes sign 
at a = ao ~ 0.65 and is negative at larger a (smaller 
fields). The implication of this result is that C(H) 
changes sign at a finite field even when the gap anisotropy 
is infinitesimally small. We analyzed the evolution of ao 
with increasing b and found [Fig. |31[b)] that oq remains 
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FIG. 4: The angular dependence of the density of states, 
Eq. (3), for 6 > 1. The form of N{cj},H) is clearly differ- 
ent from cos4(/). Additional minima or maxima correspond 
to the nodes at cos2i/) = 1/6. Panel (a) - large fields, a < 1, 
panel (b) - small fields, a > 1. The sign of the oscillating 
component still depends on whether or not a exceeds some 
6— dependent critical value. 



finite and of order one for all 6 < 1. Observe also that 
a?F{a) is of order 10~^ for all a except for the smallest 
one, i.e for 6^1, the oscillating part of C(H) is of order 
10-2 of the total C(H). 

In Fig. |3]we present A^e(H) for 6 > 1, when the gap 
along electron FSs has accidental nodes at cos 2?/; = 1 /h. 
We clearly see that the angular dependence is different 
from cos Acj) - there appear additional maxima or minima 
in iVe(H) associated with zeros of Ag (minima of |Ae|). 
These deviations from cos 40 form for b > 1 have been 
reported before 

To estimate the value of Hq in a = Hq/H, we use 
A w I8OA/VB, where B is the value of a field in 
Tesla, and take vp ~ 3.3 x lO^m/s, averaged between 
two electron bands'', and Ag ~ AmeV——. We obtain 
i?o = (0.89/c2)T. Hanaguru et aU^ extracted a smaller 
A ^ 1.7meV from their STM data. This will lead to even 
smaller Hq ~ (O.lG/c^). The value of c is not known but 
should generally be of order one. For c < 1, Hq is well 



below 9T at which the experiments are performed. In 
other words, H = 9T is deep inside the range oi H > Hq, 
when the oscillating part of C (H) has cos Acj) form with 
the maxima at — 0,7r/2, etc. Using Hq ^ IT for def- 
initeness and collecting the contributions to A^(H) from 
two hole and two electron bands, we found that the am- 
plitude of the oscillating part of C(H) is 0.028(26)2 of 
the total C(H), which for 2b = 0(1) is quite consistent 
with one percent effect observed in the experiment 

To conclude, we considered analytically, within BPT 
approximation, the variation of the field-induced com- 
ponent of the specific heat with the direction of the ap- 
plied field. We demonstrated that this scenario yields the 
cos 40 variation with the maxima at 7r/4, 37r/4, etc if the 
applied field is smaller than Hq < IT, and cos 40 vari- 
ation with the maxima at = 0,7r/2, etc if the applied 
field is larger than Hg. Both results are valid provided 
that the gaps along electron FSs have cos 2^ component, 
but no nodes. We argued that the cos 40 form of oscil- 
lating part, the positions of the maxima, and the relative 
magnitude of the oscillating component of C (H) are con- 
sistent with the experiments by Zheng et al^ performed in 
the 9T field, well above Hq. We therefore argue that the 
data on the angular dependence of field-induced C(H) 
in FeSeo.4Teo.6 are actually consistent with no-nodal ex- 
tended s— wave gap in this material. The same no-nodal 
extended s-wave gap has been extracted from STMi^, 
Andreev refiection^S, and optical conductivitj«2i data on 
FeSeo.45Teo.55- 
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